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Abstract 

The impact of large-scale magnetic fields on the interplay between the ordinary and inte- 
grated Sachs- Wolfe effects is investigated in the presence of a fluctuating dark energy com- 
ponent. The modified initial conditions of the Einstein-Boltzmann hierarchy allow for the 
simultaneous inclusion of dark energy perturbations and of large-scale magnetic fields. The 
temperature and polarization angular power spectra are compared with the results obtained 
in the magnetized version of the (minimal) concordance model. Purported compensation 
effects arising at large scales are specifically investigated. The fluctuating dark energy com- 
ponent modifies, in a computable manner, the shapes of the 1- and 2-cr contours in the 
parameter space of the magnetized background. The allowed spectral indices and magnetic 
field intensities turn out to be slightly larger than those determined in the framework of the 
magnetized concordance model where the dark energy fluctuations are absent. 



1 Motivations and goals 



Defining z ICC = 1090.51 ± 0.95 as the typical redshift of recombination according to the 
WMAP data alone (see [H [2] for the 3yr data release and 0, HI El El [7] for the 5yr data 
release) the temperature anisotropics for multipoles £ < ^z rec are practically unaffected by 
the thickness of the visibility function. The resulting angular power of the temperature 
inhomogeneities is determined by the ordinary Sachs- Wolfe effect (SW in what follows), by 
the integrated Sachs- Wolfe effect (ISW effect in what follows) and by their mutual correlation. 
The origin of the SW and of the ISW contributions is physically distinct. When the opacity 
drops suddenly, the SW term is determined by the density contrast of the photons and by the 
curvature perturbations around z rec . Conversely the ISW effect depends upon an integral 
over the conformal time coordinate r extending between r rec (corresponding to z rec ) and 
To (i.e. the present time when z — 0). The integration path determining the ISW passes 
through z\, i.e. the redshift at which the dark energy contribution is of the same order of 
the total density of matter and the geometry starts accelerating. The explicit contribution 
of the dark energy to the SW and ISW effects (as well as to their cross-correlation) has been 
scrutinized, for the first time, in [8] (see also [9]) where the effects of the late dominance of 
a cosmological term have been taken into account in a consistent semi-analytic treatment. 

The SW contribution typically peaks for comoving wavenumbers k ~ 0.0002 Mpc" 1 while 
the ISW effects contributes between fc min = 0.001 Mpc -1 and fc max = 0.01 Mpc -1 . For com- 
parison recall that the angular power spectra of the Cosmic Microwave Background (CMB 
in what follows) fluctuations are customarily assigned at a pivot scale k p = 0.002 Mpc -1 
which corresponds, for the best fit parameters of the WMAP 5yr alone (and in the light of 
the ACDM scenario^) to £ p ~ 30. Even if both contributions are reasonably separated in 
scales, the SW and ISW effects may partially compensate. An effective compensation would 
imply a suppression of the lowest multipoles (and in particular of the quadrupole) in the 
angular power spectrum of the temperature autocorrelations. In the context of the ACDM 
scenario with purely adiabatic fluctuations the compensation is partial and uneffective for 
the suppression of the quadrupole. If non-adiabatic fluctuations in the dark-energy sector 
are consistently included, the physical situation is different and a quadrupole suppression 
cannot be excluded. This possibility has been investigated especially in the context of the 
(observed) low quadrupole of the temperature autocorrelations [101 [TT1 [12] (see also [13] for 
a contending explanations of the same class of phenomena). 

The relative contribution of the SW and ISW terms does depend upon the evolution of the 
background geometry between the recombination epoch and the present time. Furthermore 
it does also depend upon the potential presence of non-adiabatic fluctuations as well as 
upon the specific way dark energy is parametrized. In the ACDM paradigm there are, by 
definition, no fluctuations in the dark energy sector. The rationale for such a statement 

1 A stays for the dark energy component (i.e. just a cosmological constant) and CDM stands for the cold 
dark matter component. 
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stems directly from the value of the barotropic index of the dark energy (i.e. w<ie = — 1)- 
The natural extension of the ACDM paradigm to the case when the dark energy fluctuations 
are dynamical is represented by the u>CDM model where w stands for the barotropic index 
of the dark energy component q In the wCDM scenario the dark energy fluctuations affect 
virtually all CMB observables and, more specifically, the best fit parameters will differ from 
the ones determined, for instance, within the standard ACDM scenario. 

A legitimate question concerns, in this context, the role of the large-scale magnetic fields. 
Indeed the ACDM paradigm served as the simplest framework for the systematic scrutiny 
of magnetized CMB anisotropics (see, e.g., pjj [15] and [HH [TTj, [18] for more recent devel- 
opments). The purpose of the present investigation is to extend the program of [HJ [T5] by 
relaxing the hypotheses of [TBI El US] on the dark energy component. What happens if the 
dominant energy density of the background is not parametrized in terms of a cosmological 
constant and, simultaneously, a magnetized background is present? For large angular separa- 
tions the interplay between the SW and the ISW effect answers satisfactorily to the previous 
question, but what happens for higher multipoles? What are the changes in the parameter 
space of the magnetized background induced by a fluctuating dark energy component? 

While there are no doubts that large-scale magnetic fields exist today in nearly all grav- 
itationally bound systems there has been mounting evidence from diverse observations in 
galaxies [TH [20] , clusters [ZD, E2] , superclusters [23] , high-redshift quasars [2H EH] that mag- 
netic fields could also have pretty large correlation scales (see [22] for a review on the subject). 
The evolution of electromagnetic fields in a plasma is subject to daily test in terrestrial lab- 
oratories and in astrophysical observations (see, e.g. [271 12H] for historic monographs on the 
subject). When plasma physics is used to interpret or even justify some astrophysical obser- 
vations the governing equations are exactly the ones used in terrestrial experiments. This 
approximation is reasonable for sufficiently small redshift. However, as we move towards 
higher redshift, the evolution of the space-time curvature cannot be neglected anymore. The 
spirit of the approach summarized in [161 El HE] is to translate the plasma dynamics in flat 
space-time to a curved background which is the one dictated, for simplicity, by the ACDM 
paradigm and by its neighboring extensions. Specific attention must be paid, in this respect, 
to the accurate treatment of the relativistic fluctuations of the geometry and to their initial 
conditions [HI [15] . 

The content of the forthcoming sections can be summarized, in short, as follows. In 
section [2] the governing equations of the system will be introduced. In section [3] the ordinary 
and integrated SW contributions will be analyzed. Section 0] deals with the initial conditions 
for the magnetized Einstein-Boltzmann hierarchy in the presence of a fluctuating dark- 
energy component. In section [5] the numerical aspects of the whole analysis have been 
collected. Section [6] contains our concluding considerations. Various technical results have 

2 To be more clear, since also other barotropic indices will intervene, we will denote the barotropic index 
of dark energy as Wdc- At the same time, since it is conventional to talk about the wCDM paradigm we will 
adhere to this convention and refrain from writing u>d c CDM. 
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been reported in the appendix to avoid excessive digressions from the main logical line of 
each section. 

2 Governing equations 

The governing equations of the pre-decoupling plasma can be written in terms of the global 
magnetohydrodynamical variables, i.e. the total current, the charge density and the baryon 
velocity (see, e.g. [29] ). In most of the astrophysical applications the space-time curva- 
ture hardly plays any role in the evolution of the plasma. Prior to decoupling, as already 
emphasized in section [TJ the space-time metric as well its relativistic inhomogeneities are 
two necessary ingredients for gauging the effects of large-scale magnetism on the CMB ob- 
servables. The reduction from the two-fluid system to the one-fluid system is discussed in 
[j~6| [T71 [18] (see, in particular, [29]) with special attention to the effects of the space-time 
curvature. The magnetohydrodynamical reduction holds, of course under various approx- 
imations which are not that dissimilar from the ones usually employed in the case of flat 
space-time. There are, however, also notable differences which can be summarized as follows: 

• since electrons and ions are non-relativistic (i.e. the plasma is cold) the system is 
not invariant under Weyl rescaling of the space-time metric and this implies that the 
corresponding evolution equations will be qualitatively different from the case of a 
non-expanding space-time; 

• even if the charge carriers are non-relativistic, the fluctuations of the geometry cannot 
be treated in the standard Newtonian approximation since the relevant modes of the 
gravitational field have wavelengths which are larger than the Hubble radius before 
matter-radiation equality §. 

With these specifications in mind, the relevant evolution equations will be written prior to 
decoupling when the electron-photon and electron-ion coupling is still strong. 

2.1 Background variables 

In a spatially flat background metric of the type g^ij) = a 2 {T)rj^ u (where r]^ u is the 
Minkowski metric) the Friedmann-Lemaitre equations take the form 

H 2 = ^^p t , H 2 -H> = 4nGa 2 ( Pt + p t ), H=-, (2.1) 
S a 

where the prime denotes a derivation with respect to the conformal time coordinate r. The 
explicit form of the energy density and of the enthalpy density appearing in Eq. ( 12. ip is 

3 The numerical value of the redshift of matter-radiation equality, denoted by z eq is given by z oq + 1 = 
3228.91[fc§n M o/(0.134)]. 
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given by: 

Pt = Pc + Pi + Pc + P 7 + Pv + Pde, (2.2) 
4 

Pt + Pt = ~{pv + P 7 ) + Pc + Pc + Pi + (Wdc + IjPde- (2-3) 

In the one-fluid limit the electron and ion energy densities form a single physical entity, i.e. 
the baryonic matter density pb = p c + Pi where p e = m e n and p x = m;n are, respectively, 
the electron and ion matter densities. The comoving concentrations of electrons and ions 
(i.e. no = a n ) coincide because of the electric neutrality of the plasma 0. The remaining 
energy densities in Eqs. (12.21) and (I2.3P parametrize the contributions of neutrinos (with 
subscript u), of photons (with subscript 7), of cold dark matter particles (with subscript c); 
finally pde denotes the dark energy density while u>d e , as already mentioned, is the barotropic 
index of dark energy. 



2.2 Fluctuations of the geometry 

Diverse gauges lead to slightly different (but mathematically equivalent) physical pictures 
of the effects discussed in the present investigation. To avoid a pedantic presentation, the 
relevant equations will be swiftly introduced in the conformally Newtonian gauge with the 
proviso that the very same equations will be expressed, when needed, in the synchronous 
frame. The recipe to move between these two gauges will now be given together with the 
appropriate definitions of the perturbed line elements. In the conformally Newtonian gauge 
the perturbed entries of the metric are given by: 

5^g 00 (k,r) = 2a 2 (r)0(fc,r), S^ 9lJ (k, r) = 2a 2 (r)^(fc, r)% (2.4) 

while in the synchronous gauge the perturbed metric can be written as: 

5i 



Si s) g l3 (k,r) =a 



2 <t) 



kikj h(k, r) + 6f r) ( hkj - 



(2.5) 



where k{ = ki/\k\. The parametrizations of Eqs. (12.41) and (12.51) are related by the appropri- 
ate coordinate transformations, i.e. 

j>(k,T) = -Z(k,T) + ^[h(k,T)+6Z(k,T)]', 

= -^{[M^) + ^r)]" + W[^r) + 6^r)]'}. (2.6) 

Consider, for instance, the curvature perturbations on comoving orthogonal hypersurfaces, 
i.e. lZ(k,r), within the description provided by the conformally Newtonian gauge: 



4 The plasma is globally neutral n\ = n c = no and the common value of the electron and ion con- 
centrations can be expressed as no — rj^nj where n 7 is the comoving concentration of photons, r/h — 
6.219 x 10~ 10 [/iofV(0.02773)][T 7 o/(2.725K)r 3 and fi b is, as usual, the critical fraction of baryonic matter. 
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where the arrow simply labels the coordinate transformation from the conformally Newtonian 
gauge to the synchronous coordinate system; the last expression in Eq. (12 .7p is obtained by 
shifting metric fluctuations according to Eq. (I2.6p . From the last equality in Eq. (12.71) it is 
also apparent that when 1Z' = 0, (,(k) = 1Z{k). In more general terms, by solving Eq. (12.71) 
in terms of £, it can be shown, after integration by part, that 

r) = TZ(k, r)-®r n) dr u (2.8) 

a(r) Jo H{ti) 

where T\ is an integration variable and where the prime denotes, as usual, a derivation with 
respect to r. Both in analytical and numerical calculations the normalization of the curvature 
perturbations is customarily expressed in terms of 71; for this reason Eqs. (12.71) and (12.81) 
turn out to be particularly useful in the explicit estimates. The same transformations of 
Eq. (12.61) can be used to gauge-transform the governing equations from the conformally 
Newtonian farme to the synchronous coordinate systeno Equation (12.81) also justifies, a 
posteriori, the perturbed form of the line element introduced in Eq. (12.51) : the factor of 6 in 
front of £ is essential if we want r) to coincide with TZ(k,r) at least in the case of the 
adiabatic mode when TZ*(k) does not depend on time prior to equality and for wavelengths 
larger than the Hubble radius. 

The present conventions as well as the whole approach slightly differ from the treatment 
of, for instance, Ref. [30] (see also [31]). In the present approach the curvature perturbations 
on comoving orthogonal hypersurfaces (i.e. TV) are simply related, in the large-scale limit, 
to the curvature perturbations on uniform density hypersurfaces (see, e.g. [32] and also 
[33j[34]). From the latter observations, most of the synchronous gauge results needed in the 
present analysis follow by making judicious use of Eq. ( 12. 8ft . The need for the synchronous 
description is also motivated since the code used for present numerical analysis is based, 
originally, on COSMICS [35] and on CMBFAST [361 [37]. The numerical code used here 
is an extension of what has been described and exploited in [T6J [T7J [18] . It is also worth 
mentioning, for completeness, that there exist approaches to cosmological perturbations 
which are fully covariant [M] and which have been also applied to the case of large-scale 
magnetic fields [39J without leading, in the latter case, to any explicit estimate or of the 
temperature and polarization angular power spectra neither in the ACDM paradigm nor in 
its neighboring extensions such as the ones analyzed in the present paper. 

Not only the metric fluctuations will change under coordinate transformations but also the inhomo- 
geneities of the sources. In particular, it can be easily shown that S^ cn ^ p = d^p — p'[(h' + 6£')/(2fc 2 )]. 
From the latter relation it also follows that a fluctuation of the energy density does not transform if its 
homogeneous background value is constant in time. 
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2.3 Evolution equations for the inhomogenities 

The Hamiltonian and momentum constraints stemming, respectively, from the (00) and (Oi) 
(perturbed) Einstein equations are given, in real space, by 

V 2 ^ - mW + if/) = A7rGa 2 [5 sPi + 5 sPdc + 5 s p B + 5 s p E ], (2.9) 



V(H(f) + itf) = —AirGa 



w E x B 



(2.10) 



where, 8 s pt = 5 s ph + 5 s p c + 8 s p v + <5 s p 7 is the density fluctuation of the fluid sources in the 
longitudinal gauge. In Eq. (12.101) the total velocity field v t obeys 

4 4 

(Pt + Pt)Vt = -puVu + gP 7 (l + Rb)v~,b + PcVc + (W de + l)pde#de; (2.11) 

as already mentioned w de is the barotropic index of the dark energy component i.e. 

Pde 2 W 'de fn 1 o\ 

w de = — , c sde = w de -— — -, 2.12) 

p de 3n{w de + 1) 

where the sound speed of dark energy has been also introduced. In Eq. (12.111) = 3pb/ (4p 7 ) 
denotes the baryon to photon ratio. Prior to photon decoupling the baryon and photon 
velocities effectively coincide with t? 7 b 

v jh ~ w 7 ~ v h , v h = ■ . (2.13) 

m e + m- s 

In Eqs. (12.91) and (I2.10p the gravitational effects of the large-scale electromagnetic fields have 
been included in terms of the comoving electric and magnetic fields E(x,t) = a 2 (r)S(x, r) 
and B(x, r) = a 2 (r)B(x, r), i.e. 

x B 2 x E<2 x $sPb x 8 s Pe / 01 a 

OsPb = - — j, dspE = - — j, o s p B = -r-, d s p E = — — , (2.14) 
87ra 4 87ra 4 3 3 



where B 2 = \B\ 2 and E 2 = \E\ 2 . The evolution equations for E and B are: 

V-5 = 0, V-£ = 4vrp q , (2.15) 
VxE + B' = 0, VxB = 4nJ + E' (2.16) 
p q + V-J = 0, p q = e(ni-n e ), J = e(n i 'u i - n e v e ). (2.17) 

The total current J has been expressed in terms of the two fluid variables, however, as in 
the case of flat-space magnetohydrodynamics, J can be related to the electromagnetic fields 
by means of the generalized Ohm law [29J : 



J = a[E + v h xB + -^ , a= y, r=-r, (2.18) 

en n e / 47r{ar ie + (4/3)[p 7 /(ra m e )jr e7 } 



6 Following exactly the same conventions established in Eqs. (|2.2|l and (|2.3|) . the various subscripts denote 
the velocities of the different fluid components. 
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where a is the conductivity; r ie and T 7e are, respectively, the electron-ion and electron- 
photon interaction rates. The three terms appearing in the Ohm's law are, besides the 
electric field, the drift term (i.e. X B), the thermoelectric term (containing the gradient of 
the electron pressur^l) and the Hall term (i.e. Jx B). It can be shown [29J that for frequencies 
much smaller than the (electron) plasma frequency and for typical length-scales much larger 
than the Debye screening length the Hall and thermoelectric terms are subleading for the 
purposes of the present analysis. The electromagnetic pressure as well as the anisotropic 
stresses enter the perturbed (ij) components of the Einstein equations: 



V>" + H{<f>' + 2V0 + (H 2 + 2U')(t> 

+^V 2 (0 - ip) = A-nGa 2 [5 s p { + 5 s p dc + 5 s p B + 5p E ] 

V 4 (0-^) 



16nGa 2 



PvV 2 <y v + p 7 V 2 a B + p 7 V 2 cr E + -pde{w d e + 1) W dc 



(2.19) 
(2.20) 



where 5 s pf = (<5 s p 7 + 5 s p u )/3, in analogy with S s pt, denotes the fluctuation of the pressure of 
the fluid components. The total anisotropic stress is given, as usual, by 



L Hv + LL iB + ii iE + ii ide) 



(2.21) 



the various subscripts denote the respective components and, in particular the electromag- 
netic contribution: 



LL iB 



1 



4vra 4 

In Eq. (12T201 the notation 

4 



B { B j - ^B 2 
3 



1 



47ra 4 



3 



4 

didjlll 3 = ^p v V 2 a v + -p 7 W B + ^p 7 W E + p de (^d e + 1) W de 



4 
3' 



(2.22) 



(2.23) 



has been adopted. The various species of the plasma either interact strongly with the plasma 
(like the elctrons, the ions and the photons) or they only feel the effects of the geometry 
(like the CDM component and the dark-energy). For large-scales the baryon-photon system 
obeys 



<Uw ^- [/xgl -v'- 



a 4 p b 



3p b 



V 2 cb + e'(9 h -9 



1)1 



5; = 4^' 



7 , 



(2.24) 
(2.25) 

(2.26) 
(2.27) 



7 In Eq. (|2.18p . following [29], what appears is the comoving electron pressure given by p c — n c T c where 
rie and T e are, respectively, the comoving concentration and the comoving temperature of the electrons. 
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where e' is the differential optical depth; furthermore, defining as Vx the velocity of the 
species X, 9x = V -Vx- The CDM and the neutrino components will evolve, respectively, as 



5' c = 3^' - 9 C , (2.28) 

e' c + w c = -v 2 </>, (2.29) 

and as 

^ = -^V 2 ^ + VV,-V 2 0, (2-30) 

K = 4^ - \e v , (2.31) 

4 3 

a' = — 6 V JF^o (2.32) 

v 15 10 " d ' v ' 

where T v z reminds of the coupling of the monopole and of the dipole to the higher multipoles 
of the neutrino phase space distribution. The latter term will be set to zero in the class of 
initial conditions discussed in the present paper but it can be relevant when magnetized 
non-adiabatic modes are consistently included. The fluctuations of the dark energy density 
obey, in the wCDM scenario, the following equation: 

S s p de + 3H(5 s p de + 5p dc ) - 3ip'(p de + p de ) + (p de + Pde) V • v de = 0, (2.33) 

where, as already mentioned, v de is the velocity of the dark energy in the conformally New- 
tonian frame whose evolution equation can be written as 

(Pde + Pde)^ e + [Pde + Pde + ^(Pdc + Pdc)]v d c + Vffyde 

+ (Pdc + Pde)V0 - (p dc + Pde)V<7 de = 0. (2.34) 

In Eq. (12.341) <7d e accounts for the possible presence of an anisotropic stress for the dark 
energy contribution. This contribution as argued in [40] (see also, in a related perspective, 
[4"Tj ) can be rather interesting in connection with dark energy models. The contribution 
of the dark energy anisotropic stress has been included, so far, to have sufficiently general 
equations. In the standard version of the uCDM scenario the dark energy has no anisotropic 
stress and therefore, for practical reasons, the anisotropic stress will be set to zero. It might 
be interesting, in the future, to relax this assumption by following, for instance, the approach 
of [10] • By introducing the density contrast for the dark energy component, Eqs. (I2.33P and 
(I2.34p can also be written as 

5' dc + 3H{c 2 sdc - w dc )5 dc + {w dc + l)0de - 3{w dc + = 0, (2.35) 

9' de + H{\ - 3c 2 sde )9 de + f deV ^ d . c + V 2 - W dc = 0, (2.36) 

{w de + 1) 

where ^de = V • v de and S s p de = p de S de . The governing equations can be easily translated in 
the synchronous coordinate system (see, in particular, Eqs. (|2.4p . (I2.5P and (I2.6p - fl2.8p ). For 
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instance the synchronous form of Eq. (12.91) can be obtained by using, directly, Eqs. (12.41) 
and (12.51) and by appreciating that 

rVf = *fp f - Pf^y^ ew = <f w - pJ-^ 1 - (2.37) 

Since pf = (p 7 + p u + p c + pb) and p dc are separately conserved (i.e. p' { = — 3H(pf + Pf ) and 
p' de = — 37i(pdo +Pde)) the result will be, as expected: 

1k 2 i - Hti = 87rGa 2 [5P Pi + 5 s (s) pde + 5 s (S) Pb + 5 & (S) Pe]- (2.38) 

In the same way all the other governing equations reported in this sections can be translated 
in the synchronous frame. 

2.4 Large-scale evolution 

The Hamiltonian constraint of Eq. (12.91) can be expressed as 

where 71 has been introduced in Eq. (12. 7p and ( [18J (see also [33], [31]) is a further gauge- 
invariant variable which measures the density contrast on uniform curvature hypersurfaces. 
The constraint (I2.39P can also be used as an explicit definition of ( whose evolution depends 
also upon the gravitating effects of the electromagnetic fields and of the Ohmic current [29] 

dC H5 Pnad | H(3cj-1) 5 

dr Pt + Pt 3(p t + pt) sPB 

E ■ J 5sP' E + 3(4 + l)H5 sPE 6 t 
3a 4 (p t + p t ) 3(p t + p t ) 3' 1 " ' 

dJ f^j , ^ 4p 7 r e7 



or V 6riQ m e 



+ / £ + 4x5 + ^-— (2.41) 
47r V erio eno / 3m e 

In Eq. ( 12.401) the contribution of non-adiabatic pressure fluctuations has been included for 
sake of completeness even if 5p nSud will not play a specific role. One of the possible extensions 
of the present work could indeed be the inclusion of non-adiabatic modes in the dark-energy 
sector (as argued, for instance, in [10] in connection with the problem of the quadrupole 
suppression). In connection with Eqs. (12.391) . (I2.40p and (12.411) few remarks are in order: 

• since both ( and 1Z are invariant under gauge transformations, they do coincide, in 
any frame, in the limit 
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• the right hand side of Eq. (I2.40p contains three qualitatively different contributions: 
the purely electric contribution (proportional to 5 s pE and its derivative), the purely 
magnetic contribution (proportional to S s pb) and the Ohmic contribution (containing 
explicitly the Ohmic current J); 

• the factor 6 t at the right hand side of Eq. (12.40)) is of order k 2 /H and it is subleading 
for typical length-scales larger than the Hubble radius; 

• Eq. (12.421) reduces, at large scales, to Eq. (12.181) : this statement can be directly 
verified by taking into account that the rate of electron-ion interactions (i.e. r ie ) is 
always much larger than the rate of electron-photon interaction (i.e. T 7e ); prior to 
decoupling both r ie and T 7e are larger than TC. 

The physical hierarchies between the electric, magnetic and Ohmic contributions follow from 
the typical length-scales of the problem and also from the frequency range dictated by the 
magnetohydrodynamical approximation (see [29J and, in particular, the left plot of Fig.l): 

W <H6 sPE « J' f ~™*e*. (2.42) 



(Pt + Pt) * a 4 (Pt + Pt) (Pt + Pt) 
The magnetic contribution appearing at the right hand side of Eq. ( 12.401) can also be written 

as 

777z 7 r-— -R 7 ^B, Rj + Rv = 1, (2.43) 

3(1 + w t ) p t 

where R v = 0.405 and R 1 = 0.595 are, respectively, the photon and neutrino fractions 
present in the radiation plasma and 

n B {x,T) = -— . (2.44) 

Recalling the parametrization of the anisotropic stress <tb, the divergence of the magnetohy- 
drodynamical Lorentz force, i.e. V • [ J x B] can be expressed as a combination of V 2 ctb and 
V 2 fi B . These vector identities are important in various analytic estimates (see [T6 | [T71 [T8] 
and discussions therein). In the case where the dark energy is simply given by the cosmo- 
logical constant the total barotropic index and the total sound speed can be written as: 

_ Pt 1 ~ 3 % 2_Pt_ a dw t _ 4 

Pt 3(l + a + %) Pt 3(1 + w t ) da 3(3o; + 4) 

where a = a/a eq and 

aA = ^ = 2246 .81 C^m) i,3 (3^\ ao = 3195.18 f^V (2.46) 
a eq V0.1326/ V0.3835/' V0.1326/ v ' 

Using Eqs. (12.451) into Eq. (12.431) we also have, in Fourier space, 

^-f ^eW = --^- 2 R^{k). (2.47) 
3(1 +w t ) p t (3a + 4) 2 



11 



3 Large-scale compensations 

The evolution of the brightness perturbations can be written as B 



A[ + {iky, + e')Ai = ip' - iky<p + e' 



Aio + fiv h + 



3/x 2 



4 



-S P (k,r) 



(3.1) 



where Sp can be expressed as the sum of the quadrupole of the intensity, of the monopole 
of the polarization and of the quadrupole of the polarization, i.e., respectively, Sp(k, r) = 
(A12 + Ap + Ap 2 ). The multipole expansion of the brightness perturbation reads, within 
the present conventions, 



Ai(fe, (i, r ) = £(-i)'(2* + 1) A«(fc, T )P t (ji), 



(3.2) 



where Pi(fi) are the standard Legendre polynomials. The line of sight solution of Eq. (13.11) 
can be written 



Ai(k, fi,T ) 



K{t) 



A 
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(1 - 3fi 2 ) _ 



-i/ix(r) 



+ 



c/re- £(r ' To) (0' + ^e-^^dr, K{t) = e'e- e(T ' T °\ 



(3.3) 



where the term —ikfM(j) has been integrated by parts^l and where JC(r) denotes the visibility 
function whose explicit form can be approximated by a Gaussian profile [4~2| H3j HU US] with 
different methods. At large scales the Gaussian can be considered effectively a Dirac delta 
function approximately centered at recombination. 

For typical multipoles i < ^z rec the finite width of the visibility function is immaterial. 
This means that for sufficiently small I everything goes as if the opacity suddenly drops at 
recombination. This implies that the visibility function presents a sharp (i.e. infinitely thin 
peak at the recombination time). Thus, since /C(r) is proportional to a Dirac delta function 
and e~ e ( T ' T °^ is proportional to an Heaviside theta function. Under the latter approximations, 
Eq. (I3.3p leads to the wanted separation between SW and ISW contributions: 



Ai(fc,/x,r ) = A[ SW) (k,fi,T ) + Aj liDVV; (/c,/i,r 



(ISW) . 



A 



(SW) 



TO 



A? sw) (fc, fi, t ) = r (0' + V/)e"^ (r) dr. 



(3.4) 
(3.5) 

(3.6) 



8 As usual, (j, = k ■ n denotes the projection of the Fourier mode on the direction of propagation of the 
CMB photon. 

9 Of course one might also wish to continue with the integrations by parts and integrate all the /i-dependent 
terms. Such a step will however produce various time derivatives of the /C(r) which are difficult to evaluate 
explicitly when the visibility function is infinitely thin. See also, in this respect, the discussion of appendix 

B 
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where it has been used that the <5 7 obeying, for instance, Eq. (I2.27p . is also related to 
Aio as 5 1 (k,r) = 4Ai (k,r) . Equations (13.51) and (13. 6p can be evaluated within three 
complementary approximation schemes: 

• [a] in the first approximation we can assume that a iec /a eq ^> 1 and that, simultaneously, 
the phase appearing in Eq. (13. 6p is r-independent (i.e. ifix(r) ~ ikfi(r — r rec )); 

• [b] in a more accurate perspective the assumption that a rec /a cq 1 can be dropped; 
after all 

«rec Z eq + 1 / h Q Q M0 \ 

"rec = = — = 3.04 — — ; (3.7) 



a eq z Tec + 1 V 0.134 

since a rec — 3 is not extremely larger then 1, significant corrections can be expected; 

• [c] finally, Eq. (13.71) can be taken into account in conjunction with the r dependence 
of the phase in the integrand of Eq. (13. 6p . 

In the absence of large-scale magnetic fields, it is sometimes practical to estimate the large- 
scale temperature autocorrelations by only keeping the ordinary SW contribution evaluated 
in the limit a rec ^> 1. The same approximation (with the same level of accuracy) can also 
be used when large-scale magnetic fields are included in the analysis. It would be incorrect 
to use different levels of accuracy for the curvature contribution and for the magnetized 
contribution. 



3.1 Simplified analysis of the compensation 

In the case labeled by [a] in the above list of items \a rec /a eq \ 3> 1 and the ordinary SW 
contribution turns out to bS 



3 R 

2i/;(k, r rcc ) - -i/j(k, tv) = 2tfj{k, r rec ) + K*(k) - ^ B {k), (3.8) 



4 

where it has been used, according to Eq. (I2.27p . that in the large-scale limit 

&y(fc, ^ec) = T rec ) + 5^k, r„) - 4^(fe, T*) + O^V). (3.9) 

Equation (13.81) then follows from Eq. (13. 9p by recalling that the initial conditions at r* are 
given, in the case of the magnetized adiabatic mode, as: 

5 7 (fc, r*) = -2Mk) - i? 7 fi B (fc), K*(k) = -ip*(k) - ^4^- (3.10) 

Since the phase appearing in the integrand of Eq. (13. 6p is estimated by positing x{j) ~ y rec 
the cancellation between the ordinary and the integrated SW contributions is maximal, and, 



3 The time labeled by is such that r <C r eq and, simultaneously, fc/7i* ~ fcr* <C 1. 
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in this sense, such an approximation can be improved: 



[2i;(k,T )-24>(k,T rcc )} 



R^n B (k) 



K*(h) + 2ij(k,r ) 



R^Q B {k) 



4 



(3.11) 



The second (intermediate) equality appearing in Eq. (13. lip has been included just to show 
the explicit cancellation at r rec . In a pure CDM model (i.e. no dark energy, no magnetic 
fields and I^mo = 1), the quantity in squared brackets in Eq. (13.1 ip can be written as: 



K,(k) + 2ij(k,r ) - 



R y fl B (k) 



K*(k) R,,n B (k) 



(3.12) 



4 54 

The estimate provided by Eq. (13.121) can be improved by taking into account the observation 
of Eq. (13.71) . In the latter case it can be shown that the SW contribution turns out to be 



SW B (a) 



5W R (« rec ) + 



R^Q B (k) 
20 



SW B (a T 



1 + 



16 



3a 3a 2 



16(vg+T - 1) 
3a 3 



12 48 32(1 - y/a + T) 
1 1 1 - 

a a z cr 



(3.13) 
(3.14) 
(3.15) 



where a rec has been introduced in Eq. (13. 7p . The result of Eq. (I3.13P corresponds to 
the approximation scheme labeled by [b] in the above list of items. The evolution of ip is 
obtained by solving explicitly the corresponding evolution equations for TZ and by recalling 
the connection between 1Z and ip (see Eq. (12.71) ). The functions reported in Eqs. (13.141) 
and (I3.15P have been derived in the appendix [A] by using, consistently, the synchronous 
gauge description. It is interesting that the synchronous description leads to an explicit 
derivation of the SW and ISW contributions which is technically different from the one of 
the conformally Newtonian gauge. By taking the limit a\ — > oo (for a fixed) in Eq. (12.451) . 
the standard CDM result is recovered. In the pure CDM case the long-wavelength limit 
implies 

ip' , H 2 -H' 



H 2 



-1> 



which can be easily solved analytically 
where 

a[a{9a + 2) -8] 



X 2 (a) 



--TZ(k,r), 

R^n B (k)l 2 (a) 
f 16[Va~TT- 1] 



(3.16) 
(3.17) 



15a 3 



3{a[a(a - 2) + 8] + 16[1 - Va + l]} 



20a 3 



(3.18) 
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In the case of the ACDM model \ip(k, r)| grows linearly as a function of the conformal 
time coordinate as a consequence of the presence of the dark-energy phase. To obtain a 
quantitative estimate Eq. (I3.16P can be used but in the case when «a and it is given as in 
Eqs. (12.451) and (12.461) . Similarly, by summing up Eqs. (12.331) and (I2.34p we shall have 



V>" + 2(W - n 2 )ip = 4xGa 2 



4 

Ode + l)6 s p de + -S s p B 



(3.19) 



In the case where w de = — 1 (which is the ACDM case), Eq. (13.191) reduces to 



1 FnrCn 2 

r = —Y—RyPR^nik). (3.20) 

In the limit flsik) — > we will have that ip{k,r) = Ci(k) + C2{k)r. The latter expression 
holds asymptotically. Direct numerical integration must be eventually employed and the 
result can be expressed as: 

iP(k,T ) = -0.4507 K*(k) + 0.1125^^0)- ( 3 - 21 ) 

Equation (I3.2ip must be compared with the value of i[)(k, r) at r rec (corresponding to z rec = 
1090.51) 

ip{k,r iec ) = -0.5882 K* (k) +0.1467 i^eO)- ( 3 - 22 ) 

The result of Eq. ( 13.22[) is also compatible with the estimates obtainable, from Eq. (12.451) . 
in the limit a a ~~ * oo: 

ip(k,T Tec ) = -0.6 7e*(fc) + 0.15 R^n B (k). (3.23) 
From Eqs. (13.211) and (13.221) we can therefore say that 

A^(Jfe) = xp(k,T ) -i){k,T IOC ) = 0.1375 72.- 0.0342 ityfiB^)- (3.24) 
Concerning the result of Eq. (I3.24p few comments are in order: 

• in the absence of large-scale magnetic fields Aif)(k) ~ 0.137 1Z*{k); the latter result 
depends upon the cosmological parameters and the quoted value refers to Eq. (I2.46P : 

• the comparative growth of the magnetized contribution is smaller than in the curvature 
case, i.e. Aip(k) ~ — 0.0342 R,yQ B (k); the sign difference stems directly from the form 
of the Hamiltonian constraint (see Eqs. (12.91) and (12.391) ) as well as from the solution 
of Eq. (I2.40p in the large-scale limit. 

It is relevant to point out that the value Aip(k) ~ 0.13772*(A;) is compatible with the results 
obtained in [ID] (i.e. 0.1472*(/c)) under the same approximations discussed here but in the 
absence of large-scale magnetic fields. 
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3.2 Transfer functions and semi- analytical estimates 



In the wCDM case the numerical integration is unavoidable even in the absence of large-scale 
magnetic fields. The semi-analytical treatment can be carried on, up to some point, and the 
obtained results are a useful complement of the numerical discussion (see section [S]) . More 
precisely the transfer function can be computed directly in terms of the evolution of the 
appropriate background quantities whose specific form will have to be studied numerically. 
Using Eq. (I3.2p the line of sight solution (I3.3P can be written as 



A K (fc, r ) 



TO 



+ 



dr 



R^VL B 
20 



Or 

SW B (a 



je[x(r)) dr 

je[x(T Iec )], 



(3.25) 



where, we recall that x{t) = k(r — r) and x(r rec ) coincides, by definition with y rec ; SW-ji(a rec ) 
and <SWB(a r ec) have been already introduced in Eq. (13.131) . The functions Tk(t) and Tb(t) 
are given by 

H(r) 



T r (t) 



K*{k) 
B(t,x) 



a 2 (r) Jo 
r a 2 (x) 



a (x) dx 



a 2 (r) 



dx + a 2 (x)dx / B(r,y)dy; 
Jo Jo 



.H{x) H(t) 

for two generic arguments x and y, the function B(x, y) is defined as 

H(x)H(y) [3ci(y) - 1] p R (y) _ 



B{x,y) 



a 2 (x) 3[w t (y) + 1] p t (y) 



(3.26) 
(3.27) 

(3.28) 



c 2 t (y) is the sound speed of the total plasma already introduced in Eq. (12.451) and here 
appearing as a function of a generic integration variable y. Equations (13.261) and (I3.27P can 
be derived, after some algebra, from Eqs. (12.71) . (12.391) . (I2.40p and (12.411) . The derivation of 
Eqs. (I3.26I) - (I3.27I) is specifically discussed in appendix IB1 

The results of Eqs. (I3.25p . (I3.26P and (I3.27P imply that the angular power spectrum for 
the temperature autocorrelations consists of three distinct terms parametrizing, respectively, 
the ordinary Sachs- Wolfe constribution, the ISW term and the cross-correlation between the 
SW and the ISW terms: 



^nf(TT) _ ^(SW) . ^(ISW) . ^(cross) 

The ordinary SW contribution turns out to be, within the present approximations 



(3.29) 



C, 



(SW) 



47T 

25 



SWl{a rec ) 



00 dk 
~k 



Vn(k)j!(kr Q ) 



+ 



7T 



100 

2tt 



00 dk 



Vn(k)j!(kr ) 



'271 f 00 dk i i 

— R 7 cos pSW B (a rcc )SW n (a vcc ) J —^V n (k)^Vn(k)j 2 (kr ), (3.30) 
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where j3 is the correlation ang 

10 

In Eq. f)3.30p we recalled that r rec <C t , implying that 
k(r — r rC c) ~ kr . The relevant power spectra and the spherical Bessel functions can be 
written ao 

/ U s n 8 -l / U \ 2(n B -l) f^T" 

' Pn(A;)= ^ B (fci:) ' ^w = v^ j m/ 2 w, (3.31) 

the integrals over the Bessel functions can be analytically performed since, as it is well known 



/ J u (ay) J lx (ay)y x dy = x , r-A r-. (3.32) 

j o 2 A r( - v+ ^ +x+i jpi |p / i^-At+A+i j 

The ISW contribution is more complicated and, as previously remarked, it cannot be com- 
puted in fully analytic terms since each term contains two time derivatives of the transfer 
functions: 



^(isw) 1R f°° dk n> dl n . fuA f T ° dT n . 

C\ J = lQn -rV n {k) —— j^kAr^d^ ——j e (kAT 2 ) 
Jo fc jo dT\ Jo ar 2 

1 Jo k Jo dT\ Jo ar 2 

+32^008^ r ^JVu(k)JVn(k) f° ^ MkArJdn H ^j e (kAr 2 ),(3.33) 
Jo k v Jo ari jo ar 2 

where Ar x = (r — ri) and Ar 2 = (r — r 2 ) and T\ and r 2 are integration variables. The last 
contribution listed in Eq. (13.291) is the cross term whose explicit expression turns out to be 

Ct° ss) = l -fsW n (a rcc ) fvAk) jf (^)MkT )UkAr)dr 
Air f°° rlk r T o /HTtd\ 

+ —R*SWl{a TCC ) J Q T V n (k) J Q I^.Jj e (kro)j e (kAr)dr 

+ ^R, cos (3SW n (a iec ) ^^(k)^Vn(k) £ {^jj^kr^kA^dr 

1 f\<jr POO Ah i i fTQ f7T!r> 

+ _Z RjCOs(3 SW B (a rec )J o —Jv^k)^Vn{k) J q ^h{kro)h{kAr)dr, (3.34) 



11 The general idea behind [HI [15] ( see a l so [El 03 HH] ) has been to include large-scale magnetic fields at all 
the stages of the Einstein-Boltzmann hierarchy and, in particular, at the level of the initial conditions. This 
means that the solutions such as the magnetized adiabatic mode and the various non-adiabatic modes are 
regular (in a technical sense) and well defined when the magnetized contribution is correctly supplemented 
by the curvature contribution. The correlation between the components is automatic already at the level 
of the initial conditions. At earlier times the correlation is also suggested, incidentally, by models where 
large-scale magnetic fields are generated during an inflationary stage [15]. 

12 The scales k p and kh are the two pivot scales at which the amplitudes of Viz(k) and of Pn(k) are 
assigned. The pivotal values of these quantities will be reminded in a moment but can be also found in Eqs. 
(3) and (5) of Ref. [16] (see also [Til [IS] for further details). 
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where kAr = k(r — r). 

Even if the time evolution of Tk(t) and 7b(t) must be obtained numerically, Eqs. (13.331) 
and (13.341) can be further simplified by anticipating the integration over k. Since the both 
Vfc(k) and Vn(k) have a power dependence upon k, the expressions appearing in Eq. (13.331) 
can be made just more explicit by taking into account that the integral (see Eq. (1C4[) ) 

AT A J v (ka) J u (kb) dk (3.35) 

is proportional to a confluent hypergeometric function [171 HEj . Using the quadratic transfor- 
mations [H], the confluent hypergeometric functions can be related to Legendre functions. 
This analysis is reported in appendix [C] (see Eqs. flC.4j) and (1C.14I) ). The expressions of 



appendix O generalize the results of [HI E] holding in the Harrison- Zeldovich limit. 

It is appropriate to remind about the conventions leading to the power spectra mentioned 
in Eq. (I3.3ip . This topic has been thoroughly discussed (see [16], (TTJ [18] and references 
therein) and will just be swiftly mentioned here. The power spectra of the magnetic fields 
are assigned exactly in the same way as the spectra of curvature perturbations: it would be 
rather weird to have a convention for the power spectrum of curvature perturbations and a 
totally different one for the magnetic power spectrum. In the wCDM paradigm (as in the 
standard ACDM case) the temperature and polarization inhomogeneities are solely sourced 
by the fluctuations of the spatial curvature 1Z whose Fourier modes obey 

(K{k)K{p)) = — V n (k)5^(k + p), v n {k)=An[ Y ) ; (3.36) 

in the ACDM case and in the light of the WMAP 5yr data alone [3 SI [5] An = (2.41 ± 
0.11) x 1CT 9 ; as already mentioned in section [1] the pivot scale k p is 0.002 Mpc -1 . In full 
analogy with Eq. (I3.36j) the ensemble average of the Fourier modes of the magnetic field are 
given by 

_ 9-7T 2 / \ * / h \ ™B-1 

(BiWBjffi) = —P tJ (k)V B (k)S^ 3 \k + p), v B (k)=A B ( K -J , (3.37) 

where Pij(k) = {k 2 5ij — kikj)/k 2 ; Ab the spectral amplitude of the magnetic field at the 
pivot scale k L = Mpc -1 [H1[T7|. In the case when n B > 1 (i.e. blue magnetic field spectra), 
A B = (2n) ns - 1 Bl/T[(n B - l)/2]; if n B < 1 (i.e. red magnetic field spectra), A B = [(1 - 
7Zb)/2](A;a/^l) "^-Bl where kx is the infra-red cut-off of the spectrum. In the case of white 
spectra (i.e. n B = 1) the two-point function is logarithmically divergent in real space and 
this is fully analog to what happens in Eq. (13.361) when n s = 1, i.e. the Harrison- Zeldovich 
(scale- invariant) spectrum. By selecting k^ 1 of the order of the Mpc scale the comoving 
field .Bl represents the (frozen-in) magnetic field intensity at the onset of the gravitational 
collapse of the protogalaxy [18] . On top of the power spectrum of the magnetic energy 
density it is necessary to compute and regularize also the power spectrum of the anisotropic 
stress and explicit discussions can be found in [T71 HE] . 
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4 Magnetized initial conditions with dark energy 



If the dark energy component is fluctuating there are, in principle various ways in which 
initial conditions of the Einstein-Boltzmann hierarchy can be set. In general terms the pre- 
decoupling plasma contains fivj^l physically different components (see, e.g. section [2]). In 
the ACDM case the dark energy does not fluctuate and, therefore, the way initial conditions 
are set does not differ from the CDM scenario insofar as the entropy fluctuations of the whole 
plasma vanish for large scales i.e. for k/Tt <C 1. The latter condition can be expressed by 
requiring that <Sy = where i and j denote a generic pair of fluids in the mixture and where 



^ = 3(0-^; 
d = -v> + 



6 



(cn) 



3( Wi + 1) 



(S) 



3{Wi + 1) 



(4.1) 
(4.2) 



are, by definition, the entropy fluctuations of the system. The second equality in Eq. (14.21) 
follows by transforming the fluctuation variables from the conformally Newtonian to the syn- 
chronous gauge. In what follows the condition of Eq. (14.11) forbids non-adiabatic fluctuations 
in the dark energy sector as well as in the fluid sector. 

If the dark energy background is consistently included in the initial conditions together 
with large-scale magnetic fields the initial conditions of the Einstein-Boltzmann hierarchy 
must be supplemented, in the conformally Newtonian gauge, by the following pair of condi- 
tions 



Sdc(k, t) 
Ode(k,r) 



'-{w de + l)^(fc) - -(w de + l)RM B (k) + O 



(4.3) 



[1 + W de )(l - 3w de ) 



k 2 



Wde 



k 2 



1 



R n B (k)-, (4.4) 

W de )(l - 3w d e) it 



where the second equation holds under the assumption that Wd c ^ —1 and that Wd e < 0. For 
the actual numerical integration it turns out to be very useful to reshuffle the dependence of 
the dark energy variables in terms of two generalized potentials defined, in Fourier space, as 

,g(k,r) (4.5) 



6de(k,T) 

S s pde(k,r) 



dyJPde 

a 
a 



1 + w dc ) 

1 + W de ) 



f( k > T ) ~ ~ w ^)d{k, t) 

f(k,T) + ~H(l-w de )g{k,T) 



(4.6) 
(4.7) 



Using Eqs. gSJ), gSP and (H2D into Eqs. (j2"^3Tl and fTXMD we do get a system for / and g: 



fj 



f + ay/Pdc 



13 This statement holds under the assumption that electrons and protons are tightly coupled by Coulomb 
scattering. If we count separately the electrons and the ions the total number of components increases to 6. 
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/' = -mf-{k 2 -- A {l-w)[2H'-{?>w ie + $)H 2 }}g 
H(l - w de ) 



(4.9) 



The corresponding equations in the synchronous gauge can be obtained just by replacing, in 
Eqs. gSD and (Oil . 

/-/, g^~g, 0^o, ^->^- (4.io) 



The resulting equations are even simpler to integrate numerically: 

g" = f, (4.ii) 

/' = -2Hf-{k 2 -^(l-w)[2H'-(3w de + 5)n 2 ]}g + a^ e ^. (4.12) 

Equations (14. lip and (I4.12p are the ones which are numerically integrated. The initial 
conditions of the remaining components of the Einstein-Boltzmann hierarchy do not change 
and are given as in [T6J [T7j. To be more specific Eq. (2.54) of Ref. [17] contains the initial 
conditions used in the numerical discussion which will be reported in section [51 



5 Numerical results 

The choice of parameters of the pivotal ACDM model in the light of the WMAP 5yr data 
alone is given by [31 HI E] 

(fi b , tt c ,tt dc , h , n s , e re ) = (0.0441, 0.214, 0.742, 0.719, 0.963, 0.087), (5.1) 

where, consistently with the established notations (see, e.g., Eqs. (13. ip and (13.31) ). e re denotes 
the optical depth to reionization. The parameters of Eq. (15.11) maximize the likelihood when 
the barotropic index of the dark energy is fixed to —1. 

The first step of the numerical analysis is to compute Tk(t) and Tb(t) when the param- 
eters are given exactly by Eq. (15. ip but the barotropic index of the dark energy assumes 
arbitrary values which can differ from the ACDM choice (i.e. w de = —1). The functions 
Tn{ T ) an d Tb(t) have been introduced in Eq. (I3.25P and control the two relevant large-scale 
contributions to the ISW effect. According to Eqs. (13.261) and ( 13.271) . Tk(t) and Tb(t) can 
be assessed once the (numerical) evolution of the background geometry has been obtained. 
Such a technique is rather useful for the analytic estimate of the asymptotic limits and it 
otherwise demands the full (numerical) solution of the background geometry; the results of 
Figs. [T] and [2] T^(r) and Tb(t) are obtained by direct numerical integration. It has been 
explicitly verified that, indeed, Eqs. (I3.26P and (13.271) are in excellent agreement with the 
numeric result once the full numerical expression of the scale factor is known. In Fig. [1] 
(plot at the left) on the vertical axis the common logarithm of the scale factor is illustrated 
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£1 =0.0441, a = 0.214,0, =0.742, h„ =0.719 £2 U =0.0441, Q =0.214, £2 J =0.742, h„ =0.71 9 

b c de b c de 




Figure 1: The evolutions of the scale factor and of Tn(y) are illustrated for different values 
of the barotropic index of dark energy Wd e - 

while on the horizontal axis the convenient time variable is given by y = 3.33 x 10~ 4 (r/Mpc) 
where r, as usual, is the conformal time coordinate. Using y as integration variable, Eq. 
( 12. II) and the related initial conditions become very simple l 14 l 

y/ufo + uj u a + u^a 1 - 31 "^, (5.2) 
—y. + V^R?/i, (5.3) 

where cur = hl(Q 7 + fl„), cu M = hl(Q h + Q c ) and u dc = ^o^de- The physical range of the 
parameters and the present normalization of the scale factor imply that ao = 1. To illustrate 
the trends induced by different values of u>d e it is useful to plot the quantities of Fig. [TJ also 
a bit outside their physical range (i.e. for a > ao). The initial condition of Eq. (15. 3p are 
imposed by making use of an exact solution of Eq. (15. 2p (as well as of the other Friedmann- 
Lemaitre equations (12.11) ) in the limit pdc — > 0. The form of the initial conditions of Eq. 
( 15. 3p implies that y^ is well before matter-radiation equality. 

Always in Fig. [TJ (plot at the right) the profile of Tn(y) is illustrated. As expected, 
around y\ (i.e. deep in the radiation dominated epoch) T-jiiyi) = 2/3. Then, as y increases, 
T-jiiy) reaches a flat plateau for intermediate times around decoupling where Tn(y) ~ 3/5 
(see Fig. [TJ plot at the right). Absent any late dominance of dark energy, T n {i)) — ► 3/5 
exactly. The (late) deviation from the CDM asymptote (i.e. 3/5) depends upon the specific 
value of Wde and determines the ISW contribution. In Fig. [2] the function Ts(y) is illustrated 

14 To avoid potential confusions, it is appropriate to remark that the integration variable of Eqs. (|5.2[) 
and (|5.3p has nothing to do with the variables x and y introduced in section [3J In Eqs. (|5.2p and (|5.3[) 
the variable y is simply given by hoy = tHq where Hq is the present value of the Hubble constant and 
ho its (dimensionless) indetermination. Conversely, in section [3] x and y have been used as (dimensionfull) 
integration variables coinciding exactly with the conformal time coordinate. 



da 
dy 
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O =0.0441,0 =0.214,0, =0.742, h =0.719 0^ = 0.0441,0 = 0.214,0, =0.742, h„ =0.719 

b c de b c de 




Figure 2: In the plots at the left and at the right the evolution of T B (y) is illustrated, 
respectively, on a linear and on a logarithmic time scale. The case labeled with Wd e = — 1 
corresponds to the best fit of the WMAP 5yr data alone in the light of the ACDM scenario. 

both on a linear time scale (plot at the left) and using a logarithmic time scale (plot at the 
right) which better represents the early evolution before matter radiation equality and across 
decoupling. The arrow appearing in the left plot of Fig. [2] marks, approximately, the present 
time, i.e. when a(yo) — Oq = 1. In short the shape of 7b (y) can be understood as follows: 

• unlike the case of Tn(y) (which reaches a constant asymptote deep in the radiation- 
dominated epoch, i.e. for y <C 1), T B (y) — ► in the limit y — > 0; 

• in the CDM paradigm Ts(y) would reach the asymptotic value 3-R 7 /20 ~ 0.089; 

• both in the ACDM and in the wCDM cases the would be asymptotes turn into an 
intermediate plateau (see Fig. El in particular plot at the right); 

• the intermediate plateau is accurately estimated by 3i? 7 /20 while the presence of a dark 
energy background is responsible for the deviation of Tb (y) from the CDM asymptote. 

The vanishing of T^(y) for y — ► is a consequence of the evolution of TZ(k,y) and, in 
particular, of Eqs. f)2.40p and (12.411) . The same observation implies that whenever lZ*(k) — > 
0, 1Z(k, y) — > for y — > 0: such an occurrence is verified (and illustrated) in the left plot of 
Fig. [3] where 1Z(k,y) is shown to vanish in the limit 1Z*(k) — > 0. 

It can be speculated that the the contribution coming from ^b(^) and TZ*(k) might 
interfere destructively leading to an overall suppression of the large-scale contribution and, 
in particular, of the quadrupole. In Fig. [3] (plot at the right) this possibility seems excluded 
unless ^b(^) — 7l*(k). But such an extreme value would jeopardize the agreement of the 
theory with the acoustic region. 
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Figure 3: The evolution of lZ(k,y) in the case 1Z*(k) — > (plot at the left). The evolution 
of ip(k, y)/lZ*(k) for different values of the magnetic energy density. 



To scrutinize this point in a model-independent perspective, consider the condition 
fi B (/c) ~ 1Z*(k) in loose terms and postulate that the amplitude of curvature perturba- 
tions matches approximately the amplitude of the power spectrum of the magnetic energy 
density Neglecting, for simplicity, the dependence upon the spectral indices we get a typical 
amplitude for the comoving amplitude of the magnetic field intensity which is of the order 
of -Bl — 22.68 nG where it has been assumed that An = 2.4 x 10 -9 , k^ = 1/Mpc and 
k p = 0.002 Mpc -1 . The putative value of about 23 nG is indeed quite large and it would 
totally disrupt the structure of the acoustic oscillations. This aspect can be appreciated 
from Fig. H] where a magnetic field of 10 nG already jeopardizes the observed features of the 
temperature autocorrelations 0. 

In the minimal scenario [16J, [52] , the parameter space of the magnetized CM B anisotropies 
can be described in terms of the magnetic spectral index n-Q and the regularized magnetic 
field amplitude B L . These two quantities and their relations with Q B and <7b (introduced 
in section [2]) have been thoroughly discussed in previous papers [HH [T71 HE]. The basic 
advantages of such a parametrization have been illustrated at the end of section [3] (see, in 
particular, Eq. f!3.37|) and Eq. (13.361) for the analog conventions employed in the definition 
of the power spectra of curvature perturbations). In Fig. H]the parameters are fixed to the 
best fit of the WMAP 5yr data alone in the light of the ACDM paradigm (corresponding to 
Wde = — 1)- Always in Fig. H]the barotropic index of the dark energy is increased from — 1 
to —0.6 while the magnetized background is switched on. The magnetic field intensity has 
been chosen to be rather large (i.e. O(10 nG)). Examples of blue (i.e. n-Q = 1.5) and red (i.e. 

15 The threshold field of about 23 nG is obtained by requiring that Va(k) ~ Vfi(k). The latter condition 
implies, from Eq. (f3~3lj) , that £ B ~ An- Since £ B oc n| L and fl B L = 9.536 x 10" 8 (B L /nG) 2 , we do get 
Sl ~ 22.6nG. This argument holds rigorously in the scale-invariant limit and it can be easily generalized to 
the case n B ^= \. 
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Figure 4: The angular power spectra of the CMB observables in the light of the tuCDM 
scenario when large-scale magnetic fields are included. In all the plots the barotropic index 
of the (fluctuating) dark energy background has been fixed to w<ie = —0.6. 

n-Q = 0.5) spectral indices are illustrated. From the TT correlations in semilogarithmic 
coordinates (top-right plot) as well as from the TE and EE correlations the general interplay 
of the dark energy and magnetized backgrounds is more clear and can be summarized as 
follows: 

• in the large-scale domain (i.e. £ < 50) both the dark energy and the magnetized 

16 The TT power spectra denote, as usual, the autocorrelations of the temperature. The TE power spectra 
are the cross-correlations between temperature and poalrization. The EE power spectra denote the polar- 
ization autocorrelations. Within the conventions employed in the present paper, the definitions of the TT, 
TE and EE power spectra in terms of the solutions of the heat transfer equations can be found in p~7] (see 
also [18] and references therein) . The B-mode polarization induced by Faraday rotation [501 EI] wn l n °t be 
specifically addressed in this context. It has been actually shown in [SU] (second paper of the list) that the 
present data on the B-mode autocorrelations only allow for upper bounds on the Faraday-induced B-mode. 
An exception to this statement are the QUAD data (53[ [54] [55] but, at the moment, its not clear to what 
extent the results are contaminated by systematics [53] . 
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contribution augment the power in the TT angular power spectra; 

• possible interference effects leading to an overall suppression of the lower multipoles 
are excluded if the acoustic oscillations are to be reproduced correctly; this statement 
holds in the absence of any non-adiabatic contribution (i.e. 5p n8i d = in Eq. (12.401) ) 
but it might change if <5p na d 7^ 0; 

• the interplay between the magnetized and the dark energy backgrounds leads to a 
combined distortion of the peaks in the TT, TE and EE angular power spectra; since 
the level of distortion increases with the multipole the higher peaks look also shifted 
in comparison with the patterns exhibited by the underlying ACDM model with the 
same parameters and in the absence of any magnetic field. 

The value Wd e = —0.6 chosen in Fig. H]is rather extreme insofar as it is excluded by current 
fits to the WMAP 5yr data (either alone or in combination with other data sets). It is 
however useful highlight some general trend induced by the dark energy fluctuations. 

Figures [TJ E] and [3] do also contain examples where Wde < — 1- When analyzing the 
cosmological data sets in the light of a fluctuating dark energy background it can happen 
that the central value of the barotropic index maximizing the likelihood for a specific data 
set gets smaller than —1. In the latter case future singularities are expected (see for instance 
[4~9]). In what follows the cases Wde < — 1 will be considered for completeness. Indeed the 
central values of the barotropic index is still debatable and does depend upon the data sets 
which are combined in the analysis. For the purposes of the present paper, the situation is 
summarized in three forthcoming paragraphs. 

By analyzing the WMAP 5yr data alone in the light of the u>CDM model we do get 
Wde = — l-06°:g 1 42 while the remaining parameters are determined to be 

(fib, f2 c ,f2do, ho, n s , e re ) = 

(0.046°_a, 0.221^ 7 82 , 0.73 3 , n , 0.74^ 5 14 , 0.963^, 0.086^ 7 16 ); (5.4) 

in the case of the parameters of Eq. (15.41) the corresponding value of the normalization of 
the curvature perturbations is still An = 2.41 x 10~ 9 at the pivot scale of k p = 0.002 Mpc -1 ; 
the value of An does not change in comparison with the corresponding value determined in 
connection with the parameters of Eq. (15.11) . 

By combining the WMAP 5yr data with the data stemming from the analysis of the 
baryon acoustic oscillations (BAO in what follows) [56J, [571 EE] the preferred range of values 
of Wde becomes — 1.15° o^; the remaining parameters turn out to be [3j HJ E]: 

(fib, fic, fide, ^-0; n s, Cre) = 

(0.0419^56, 0.213°^, 0.745™",, 0.739™"*, 0.958™^, 0.083^ 6 16 ); (5.5) 

the value of An determined in connection with the data of Eq. (15.51) increases in comparison 
with the corresponding values holding in the case of Eqs. (15. ip and (15 .4p and it is given by 
A n = 2.48 x 10~ 9 . 
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Finally, by combining the WMAP 5yr data with the BAO and with all the supernovae 
[59| [60j [61] the central value of Wd e maximizing the likelihood gets above —1 and it is 
u>de = — 0-972 a Q 6 1 60 ; the remaining parameters are given, in this case, by [3J, SI E] 

(fib; fiofide; h-Q, n s , e rc ) = 

(0.0467^18, 0.231^ u , 0.722°;^, 0.697^ 4 14 , 0.962^, 0.085 T 6 16 ); (5.6) 

again the value of An determined in connection with the data of Eq. (I5.5P diminishes 
in comparison with the corresponding values holding in the case of Eqs. (I5.5P and it is 
A n = 2.43 x KT 9 . 

The data reported in Eqs. (15. 4p . (15.51) and (15. 6p will be referred to as, respectively, 
WMAP 5yr alone, WMAP 5yr + bao, and WMAP 5yr + bao + snail. Bearing in mind the 
latter shorthand notations it is now appropriate to scrutinize how the simultaneous presence 
of large-scale magnetic fields and of the dark energy fluctuations affects the temperature and 
the polarization observables. This aspect is illustrated in three of the forthcoming figures, 
i.e. in Figs. [5], [6] and [7] where the magnetized background is studied together with the dark 
energy background for the reference models discussed in Eqs. (15.41) . (15.51) and (I5.6p . In Fig. 
E]the parameters of the w;CDM model have been chosen as in Eq. (15 .4p with w^e = —1.06 
corresponding to the best fit in the light of the WMAP 5yr data points. It is superficially 
clear that the string of parameters of Eq. (15.11) does not differ, qualitatively, from the ones 
of Eq. (15 .4p . As expected, the inclusion of the magnetized contribution distorts the higher 
peaks. Three representative values of the parameter space of the magnetized uCDM scenario 
are illustrated: 

• (i) the curves with n-Q = 1.5 and B^ = 10 nG are representative of the region with 
high magnetic field intensity and blue spectral index (i.e. tib > 1); 

• (ii) the curves with hb = 0.5 and B^ = 10 nG are representative of the region with 
high magnetic field intensity and red spectral index (i.e. hb < 1); 

• (Hi) the curves with hb = 1.5 and B^ = 5 nG are representative of the region with 
intermediate magnetic field intensities and blue spectral index. 

The regions of the parameter space corresponding to (i) and (ii) are roughly incompatible 
with the observed angular power spectra as it can be easily argued by superimposing the 
curves to the binned data of the TT and TE correlations [3j HI [5] . Furthermore, the values 
of the parameters of regions (i) and (ii) lie in a portion of the parameter space which has 
been excluded by estimating the parameters of the magnetized background in the light of the 
magnetized extension of the standard ACDM scenario. Does the inclusion of a fluctuating 
dark energy background make the difference? The answer to this question will be the subject 
of the remaining considerations of this section. Notice also that the choice of the parameters 
corresponding to (Hi) seems to be qualitatively compatible with the observations. One of 
the purposes of the following discussion will be to make such a statement more quantitative. 
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Figure 5: Large-scale magnetic fields are included together with dark energy fluctuations in 
the case of the uGDM scenario (see Eq. (j5.4p for the list of the parameters). In the TT and 
TE angular power spectra the binned points of the WMAP 5yr data are also reported. In 
all plots the fluctuating dark energy contribution has been included and, consistently with 
Eq. (JEU), w de = -1.06. 

As a last remark concerning Fig. [51 the value of the barotropic index of dark energy 
(i.e. w^e — —1.06) is rather close to the ACDM determination and this occurrence excludes 
interference effects such as the ones observed in Fig. H] where w<ie differs from —1 appreciably 
(i.e. Wde = —0.6). In Fig. [6] and [7] the dark energy background and the other hjCDM 
parameters are fixed, respectively, as in Eqs. (I5.5p and Q5.6j) . The considerations suggested 
by Fig. H]are confirmed, in some sense, also by Figs. [5]and[6j Extreme values of the magnetic 
field intensity O(10 nG) are incompatible with the observations both in the case of red and 
blue spectral indices. 

The joined two-dimensional marginalized contours for the various cosmo logical param- 
eters identified already by the analyses of the WMAP 3yr data are ellipses with an ap- 
proximate Gaussian dependence on the confidence level as they must be in the Gaussian 
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Figure 6: Large-scale magnetic fields are included for the wCDM best fit when the WMAP 
5yr data are combined with the ones stemming from the baryon acoustic oscillations. As 
in Fig. the fluctuations in the dark energy background have been taken into account; 
Wde = —1.15 (as it must be from Eq. (15.51) ). 

approximation (see, e. g. [I]). The parameter space of the magnetized CMB anisotropics 
can then be scanned in the presence of dark energy fluctuations by using the same technique 
discussed in [TBI EH] an d applied to the case of the ACDM paradigm. In Fig. [8] the filled 
(ellipsoidal) spots in both plots mark the minimal value of the \ 2 (i- e - Xmin) obtained in the 
context of the WMAP 5yr data alone (see, e.g., Eq (15.41) ); the closed curves represent the 
likelihood contours in the two parameters tib and B^. In both plots of Fig. [8] the boundaries 
of the two regions contain 68.3 % and 95.4 % of likelihood as the values for which the x 2 has 
increased, respectively, by an amount A% 2 = 2.3 and A% 2 = 6.17. In the plot at the left 
of Fig. [JJ the data points used for the analysis correspond to the measured TT correlations 
contemplating Att = 999 experimental points from I = 2 to I = 100. In the plot at the 
right the data points include, both, the TT and TE correlations and the total number of 
data points increases to Apr + Ate = 1998. The values of n-Q and B^ minimizing the \ 2 
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Figure 7: Large-scale magnetic fields are included for the w;CDM best fit when the WMAP 
5yr data are combined with the baryon acoustic oscillations and with all supernova data. 
As in Figs. and [6] the fluctuations in the dark energy background have been taken into 
account; w de = —0.972 (in agreement with Eq. (15. 6p ) 

when only the TT correlations are considered turn out to be n-Q = 1.883 and B^ = 4.982 nG; 
in this case the reduced \ 2 is 1-094. When also the TE correlations are included in the 
analysis the reduced x 2 ; i- e - Xdof? diminishes from 1.094 to 1.033 and the values of n B and 
i?L minimizing the x 2 become n-Q = 1.913 and B^ = 5.163 nG. The points of the parameter 
leading to Xmm maximize the likelihood in the Gaussian approximation (see [TB] and [T7]). 

In a frequentist approach, the boundaries of the confidence regions represent exclusion 
plots at 68.3 % and 95.4 % confidence level. In Fig. [8] (plot at the right) the inner dashed 
curve corresponds to the 68.3 % boundary while the outer dashed line corresponds to the 
95.4 % boundary as determined in the context of the magnetized ACDM scenario and with 
the very same set of data [16J. By comparing the dashed and full contour plots in Fig. [8] 
(right plot) the shapes of the excluded regions are compatible in the two cases. The overlap 
seems even to increase by increasing the confidence range. At the same time the addition of 
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n B = 1. 883, B L = 4.982 nG n B = 1.913, B L = 5.163 nG 




Figure 8: Likelihood contours in the plane (tib, -Bl)- in both plots the inner curve cor- 
responds to Ax 2 = 2.3 (i.e. 68.3% or la of likelihood content) while the outer curve 
corresponds to A% 2 = 6.17 (i.e. 95.4% or la of likelihood content). 

a fluctuating dark energy background pins down systematically larger values of the magnetic 
field parameters. More specifically, by looking at the parameters corresponding to Xmin' we 
can say that 

(n B ,S L ) AC DM = (1-598, 3.156nG) -> (n B ,B L ) wC BM = (1.883, 4.982 nG), (5.7) 
{n B , £ l )acdm = (1-616, 3.218nG) -> (n B , B L ) wCBU = (1.913, 5.163 nG), (5.8) 

where Eq. ( 15.7ft corresponds to Ntt — 999 (i.e. only the unbinned points of the TT 
correlations are used) while Eq. (15.81) does correspond to iV TT + iV TE = 1998 (i.e. the 
unbinned points of the TT and TE correlations are used). 

Equations ( 15. 7ft and (15.81) show that the magnetic field parameters minimizing the \ 2 
(and maximizing the likelihood) do (slightly) increase when the dark energy component is 
allowed to fluctuate. Of course the slight numerical difference does not change the physical 
conclusions obtained in the context of the magnetized ACDM class of models: the preferred 
region is for blue spectral indices and moderate magnetic field intensities. In Fig. [9] the 
likelihood contours are obtained when the best fit parameters are determined from the wCDM 
parameters arising when the WMAP 5yr data are combined with the data stemming from the 
baryon acoustic oscillations (see Eq. (15.51) ). In this case, again, the shapes of the likelihood 
contours are qualitatively a bit different but quantitatively compatible with the cases of Fig. 
|8] and of Eq. [91 As expected a decrease of lowers the value of the regularized magnetic 
field intensity B^. This effect is illustrated in Fig. [10J In the right plot, the likelihood 
contours are computed when /cl = 0.1 Mpc -1 ; the parameters of the u>CDM paradigm are 
determined from Eq. (|5.4|) . In the plot at the left the same analysis is performed for the 
ACDM case when, again, the other parameters are determined from the WMAP 5yr data 
alone. The exclusion plots of Figs. 151 191 and [TOl explain, a posteriori, why in the plots of Fig. 



30 



n B = 1.674 , B L = 3.377 nG n B = 1.697 , B L = 3.461 nG 




Figure 9: Same quantities as in Fig. [8] but in the case when uCDM scenario is analyzed in 
the light of the WMAP 5yr +bao data (see Eq. fl53|) V 
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Figure 10: The effect of change in the magnetic pivot scale kj, is illustrated in the case of 
the WMAP 5yr data alone in the ACDM case (left plot) and in the uCDM case (plot at the 
right). 

El [6] and [7] the minimal amount of distortions was provided by the dot-dashed curves. The 
dot-dashed curves corresponded, in those figures, to the region (Hi) i.e. the region where the 
moderate magnetic field intensities corresponded to blue magnetic spectral indices. 

6 Concluding remarks 

The analysis of magnetized CMB anisotropies has been extended, for the first time, to the 
case when the dark energy background fluctuates. This class of models, conventionally 
labeled as magnetized wCDM, generalizes the magnetized ACDM scenario. The ordinary 
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and integrated Sachs- Wolfe effects have been analyzed in the light of the wCDM paradigm. 
In the minimal realization of the scenario (i.e. in the absence of non-adiabatic components 
in the initial conditions) the contribution of the magnetized background is not sufficient to 
compensate the effect induced on the adiabatic mode by the late-time dynamics of the dark 
energy background. Possible compensating effects arise for excessively large values of the 
magnetic field intensity; these values are already constrained by the structure of the acoustic 
oscillations. 

Since the modified initial conditions of the Einstein-Boltzmann hierarchy allow for the 
simultaneous presence of dark energy perturbations and of large-scale magnetic fields, the 
morphology of the distortions arising in the magnetized wCDM scenario has been compared 
with what happens when the dark energy background is not dynamical. The modifications 
induced by the fluctuating dark energy component on the shapes of the 1- and 2-a contours in 
the parameter space of the magnetized background have been computed. At 95 % confidence 
level, the allowed spectral indices and magnetic field intensities turn out to be systematically 
larger than those determined in the framework of the magnetized ACDM paradigm. 
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A Synchronous frame 

It is understood that all the quantities appearing in this section of the appendix are evaluated 
in the synchronous frame. The evolution equation of the intensity brightness perturbations 
are given by 



A[ + (ik/i + e')A : = — ^-(h + 60' + i 
v h + Hv h + — (3zA n + v b ) + ik — — . 



Aio + nv h + 



[l-3fi 2 )S P (k,r) 



(A.l) 
(A.2) 



The formal solution of Eq. (IA.1I) can be written, after integrating once by parts, as: 
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where, as in the bulk of the paper, x(r) = A;(r r). The visibility function is defined in terms 
of the optical depth are defined as in Eq. (13. 3p . The integration by parts can be further 
used in Eq. (EQ) : 
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(A.4) 



(A.5) 



but in spite of the simplifications obtained in Eqs. ( 1A.4I) and (1A.5I) . Eq. (1A.3I) is still the 
most direct for separating the ordinary and the integrated Sachs- Wolfe contribution in the 
presence of a sudden drop in the opacity at r rec . In this case Eqs. (1A.3I) implies 

Ai(k,fi,T ) = Af W) (fc,/i,r ) + Aj ISW) (fc,/i,r ), 
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Recalling that, by definition, y T . 
rearranged as 

Af W W,r ) = 
where the evolution equation 
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has been used explicitly. In Eq. flA.lOp the various anisotropic stresses have been neglected; 
Eq. (1A.10I) is the synchronous version of Eq. (I2.20p . as it can be easily verified by using 
the procedure swiftly illustrated in section [TJ Recalling that Ai = 4<5 7 , the solution of 
synchronous gauge version of Eq. (12.271) implies 



A 



... ■-— + £ + 0(*V). 



(A.ll) 



To obtain the explicit form of Eq. (|A.9[) we need the evolution of r) and h(k, r) across 
equality. By using Eq. (12.81) and (12.381) we obtain 



ak,T) = n*(k)- R ^ {k) g 1 (a) 



(A.12) 
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where a = a/a eq = [(t/tx) 2 + 2(t/ti)] solves Eq. (12.11) across equality and where the five 
functions Qi(a) are given by: 
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Using Eqs. flA.15p - flA.19p the ordinary SW contribution can be written as 
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where SWn{oirec) and 5WB(« r ec) are the two functions appearing, respectively, in Eqs. (13.141) 
and fl3TT5]) . 



B Explicit derivation of Tn(r) and Tq(t) 



Consider Eq. (12. 7p in terms of the conformally Newtonian variables. The first equality can 
be written, for immediate convenience, as 
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where we only passed from the conformal time coordinate r to the cosmic time coordinate 
t and we recalled that dt = a(r)dT; it is also practical to use the well known identities 
H = (H' — H 2 )/a 2 and H = H/a. From Eq. (IB.lj) by integrating once with respect to t the 
solution for ij)(k,t) then becomes 

M,t) = -K(k,t) + — / -rrdt' + - K(k,t')a(t') dt' , (B.2) 

a Jo H a Jo 

where we integrated by parts and used that 

From Eqs. (I2.39j) . (12.401) and (12.411) . neglecting the electric and Ohmic contributions, 7l(k, t) 
is given by: 

Tl^-^P^+fcl-l)^ (B.4) 
whose formal solution can be written as: 



W io Pt + Pt Pt + Pt V st 37 v ; 



H(t ')5 PoaA (k,t') df + /•*' H(t")S eP B(k , t") ( 2 1 
Pt + Pt ft + Pt 

If we now insert Eq. (IB .5 1) into Eq. (IB. 21) and go back to the conformal time coordinate we 
can write 

ip(k,r) = -n,(k)T n (r) + i? 7 n B (A:)T B (r) (B.6) 

where T^(r) and 7b (r) are given by Eqs. (13.261) and (13.271) . For sake of simplicity and for 
comparison with other treatments consider, for instance, the case of the ACDM scenario 
where 

r /3 M 2/3 
a(t) = ax sinh f -H t J . (B.7) 

In the cosmic time coordinate T-^(t) can be written as 



T n {t) = 1 - ^ f a(t')dt' (B.8) 



a(t) Jo 

It is now easy to show that, for t — > 0, a{t) ~ t 2 / 3 and Tn(t) ~ —(3/5). For t — > oo 
Tn(t) ~ e~ H ° l . Unfortunately, in the most general situation of the u>CDM scenario, the 
semi-anaytic estimates are not that simple and the numerical evaluation of Th(t) and 7b (r) 
is often mandatory (see, for instance, the initial part of section [5]). 



C Generalization of certain integral representations 

In the present study there is often the need of generalizing certain integral representations 
of the power spectrum of the ISW and of its correlation with the SW effect. For instance, in 
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the simplest case of adiabatic fluctuations we shall have that the autocorrelation if the ISW 
and its cross correlation with the SW term boils down to the following pair of expressions: 



a 



(ISW) 



where 



a 



fe(a,b, n s 



(cross) 



1 

5 

8n 2 A n k 



T, c dr\ 

to dTn 

roc. dl~ 



dJn 

r„c dTo 



dT 2 f e {T 1 ,T 2 ,n s 



dr 



TO 



1— n s 



'a b 



k ns 3 J e+ i/ 2 {ka)J e+1/2 (kb). 



(C.l) 
(C.2) 

(C.3) 



The expressions of Eqs. fIC.ip and (10. 2[) are just illustrative since the same problem arises 
in the case of the transfer functions of the magnetic component. To evaluate expressions 
like the ones of Eqs. (IC.ip and (1C.2|) one might choose to perform first the (numerical) 
integrations over the conformal time coordinate and then the integrations over the modulus 
of the wavevector k. The idea is to reverse this order and integrate first over k. In doing so 
we are led, in general terms, to integrals like 

k- x J u (ka)J u (kb)dk = 
a v r(u + 



1-A 
2 



2 X b u-X+l T ( A±l )r(i/ + l) 



1-A 1-A 



1; b 2 . 



(C.4) 



+ 1/2) and A = 3 



In terms of the illustrative examples of Eqs. (1C.1I) and (1C.2|) . v 
Equation (1C.4I) can be found, for instance, in Ref. [47] (see formula 6.574, page 675). 

It is more practical, for the purpose of approximate evaluations, to relate the hypergeo- 
metric functions to the associated Legendre functions. This was the strategy followed in the 
first attempt to evaluate analytically the ISW contribution in the ACDM model [8]. The 
latter analysis has been performed is the case of (exactly) scale-invariant contribution. In 
spite of the fact that the explicit dependence upon the spectral index is unimportant for 
orders of magnitude estimates, the habit of computing analytically the ISW contribution 
for the exactly scale- invariant case persists [TTJ. In what follows the integral representation 
of the two-point function for the ISW will be directly related to the associated Legendre 
functions but for a generic spectral index. 

Consider an hypergeometric function parametrized as F(a, /3; 7; z). It can be shown that 
if two of the numbers (7 — 1), (a — (3) and (a + (3 — 7) are equal (or if one of them is equal to 
1/2) then there exist a quadratic transformation which allows to connect the hypergeometric 
function to a Legendre function whose argument will be suitable for (approximate) semi- 
analytic integrations. Indeed in the case of Eq. flC.4j) a — v + (1 — A)/2, /3 — (1 — A)/2 and 
7 = u+1; this observation implies that (a—0) and (7 — 1) are both equal to v. Consequently, 
the general quadratic relation [48J 

1 + y/T= 



F(a, (3; 2/3; z) 



-2a 



a, a - (3 + -; (3 + -; 



1 - yi^ 

1 + VT^z~ 



(C.5) 
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implies, in the case of Eq. (IC.4I) 



1-A 1-A 



2 ' 2 ' ' & 2 
It also follows that tha10 

/ k~ x J u (ka)J u (kb)dk 
Jo 



b + a 



2v+l -A 



iH — , 2^ + -; 2i/+l; 



4a6 



(a + 6) 2 



(C.6) 



r i/ + 



1-A 



2 A (a + 6) 



2i/+l-A 



r(i^)r(^ + i) 



1-A 1 4ab 



(C.7) 



The hypergeometric function appearing in Eq. (1C.7j) can be further transformed according 
to the general relation (see |48j , p. 560) 



F(a, 2/3; z) 

In our case Eq. (1C8I) implies 

1-A 



a a + 1 



2' (2-z) 



(C.8) 



^ H , v + -; 2v + 1; 

2 2 



4a6 



a 2 + b 2 - 
{a + b) 2 _ 

Using Eq. (1C9|) we also have that: 



v 1-A i/ A 3 4a 2 6 2 
2 + 2 ~ 4 + 4 ;Z/ + 1; (a 2 + 6 2 ) 



(C.9) 



£T A J v (ka) J u (kb)dk 



a v b v 



1-A 
2 



2A (« 2 + ^ +(1 - A)/2 r(A±i)r(, + i) 



x F 



2 ^ 1-A 4a 2 6 2 
2 " 4 + 4' 2 + ^F ;Z/+ ' (a 2 + 6 2 ) 2 



(CIO) 



the hypergeometric function appearing in Eq. (IC.lOj) can be finally related to the associated 
Legendre functions according to the general relation 



n + m 1 m + n 



3 1 



1+ 2 '2 + 2 ' n+ 2'^ 2 



r m+n+l nn+1 



(z 2 - l) m / 2 v^F T(m + n + 1) 



(an) 



17 We correct here a typo appearing in the formula 6.576 (number 2) page 576 of Ref. [47]. In the numerator 
of the formula we read 2 V which should be instead b u . Equation (|C.7j) represents the correct version of the 
quoted formula where b v (and not 2") correctly appears in the numerator. 
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where Q™(z) are the associated Legendre functions of second order. Matching Eq. (IC.11I) 
with the expression appearing in Eq. (1C.10I) we shall have 

v A 3 v 1-A , 4a 2 b 2 
--T + T, - + ^— ;^ + l; 



2 4 4' 2 4 ' ' (a 2 + 6 2 ) 2 . 
I 2a6 J W + b 2 ) r ^ + i±A^ ^ 1/2 V 2ab Y [ ^ LZ) 



where we used the fact that 



2 



Q; m (y) = e- 2 — (C13) 
1 (n + m + 1) 



Thus, in conclusion, we shall have that 



OC 



k- x J u (ka)J u (kb)dk 

l-a\ 

~~n I „-i7rA/2 

-,A/2 



( a 2_62)A/2 r^ + V j e- OTA / 2 ^ A/2 fa 2 + b 2 
t(^V(u + 



2VS rfnArL^aA Qu - l/2 \ 2ab J' (ai4) 



2/1 2 



This expression generalizes, to the best of our knowledge, the analog expressions derived 
in [5] in the case of scale-invariant curvature perturbations (corresponding, in the present 
parametrization, to A = 3 — n s , i.e. A = 2 for n s — 1 and v = i + 1/2). The very same 
expression can be used to integrate the magnetized contributions to the ISW. In the latter 
case, unlike the expressions of Eqs. (IC.lj) and (1C.2|) . the various terms will contain (at least) 



one Tb(t). For instance, in the case of the magnetized contribution to the ISW there will be 
two derivatives of Tb(t) and the new A will be given as A = 4 — 2ne- All the different terms 
arising in Eqs. (13.331) and (13.341) can be explicitly integrated over k by using the techniques 
reported in this part of the appendix. 
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